Dark Matter annihilations in Pop III stars 
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We study the impact of the capture and annihilation of Weakly Interacting Massive Particles 
(WIMPs) on the evolution of Pop III stars. With a suitable modification of the Geneva stellar 
evolution code, we study the evolution of 20 and 200 Mq stars in Dark Matter haloes with den- 
sities between 10 s and 10 11 GeV/cm 3 during the core H-burning phase, and, for selected cases, 
until the end of the core He-burning phase. We find that for WIMP densities higher than 5.3 



io lu K"/io- 



GeV cm the core H-burning lifetime of 20M© and 2OOM0 stars exceeds 



the age of the Universe, and stars are sustained only by WIMP annihilations. We determine the 
observational properties of these 'frozen' objects and show that they can be searched for in the 
local Universe thanks to their anomalous mass-radius relation, which should allow unambiguous 
discrimination from normal stars. 



In the Standard Cosmological Model, the matter den- 
sity of the Universe is dominated by an unknown com- 
ponent, approximately 5 times more abundant than 
baryons, dubbed Dark Matter (DM). Among the many 
DM candidates proposed in the literature, Weakly In- 
teracting Massive Particles (WIMPs), i.e. particles with 
mass O(100) GeV and weak interactions, appear particu- 
larly promising, also in view of their possible connection 
with well motivated extensions of the Standard Model of 
particle physics (see Ref. [l| for recent reviews on particle 
DM, including a discussion of ongoing direct, indirect and 
accelerator searches). Despite their weak interactions, 
WIMPs can lead to macroscopic effects in astrophysi- 
cal objects, provided that they have a sizeable scattering 
cross section off baryons. In this case, in fact, DM par- 
ticles traveling through stars can be captured, and sink 
at the center of the stars. Direct searches and astrophys- 
ical arguments, however, severely constrain the strength 
of DM-baryons interactions (see e.g. Ref. Q and refer- 
ences therein). Since the capture rate is proportional to 
the product of the scattering cross section times the lo- 
cal DM density, large effects are thus expected in regions 
where the DM density is extremely high (this was already 
noticed in the context of the so called 'cosmions' 0]). Re- 
cent progress in our understanding of the formation and 
structure of DM halos has prompted a renewed interest 
in the consequences of DM capture in stars, in particular 
in the case of White Dwarfs fj], compact objects Q and 
main sequence stars @ at the Galactic center, where the 
DM density could be extremely high 

Alternatively, one may focus on the first stars, which 
are thought to form from gas collapsing at the center 
of 10 6 - 10 8 Af o DM halos at redshift z < 10 - 30. In 
fact, Spolyar, Freese and Gondolo @ have shown that the 
energy released by WIMPs annihilations in these mini- 
halos, during the formation of a proto-star (thus even 
before DM capture becomes efficient), may exceed any 
cooling mechanism, thus inhibiting or delaying stellar 
evolution (see also Ref. [§]). The formation of proto-stars 
with masses between 6M Q and 600M Q in DM halos of 



10 6 Mq at z=20, can actually be delayed by ~ 10 3 - 10 4 
yrs Once the star forms, the scattering of WIMPs 
off the stellar nuclei becomes more efficient and a large 
number of WIMPs can be trapped inside the gravita- 
tional potential well of the star. The WIMPs luminosity 
can overwhelm that from nuclear reactions and therefore 
strongly modify the star evolution (Tf], [T§] , and the core 
H-burning phase of Pop III stars, in DM halos of density 
of 10 11 GeV cm -3 , is substantially prolonged, especially 
for small mass stars (M* < 4OM ) [10| . 

In this letter, we perform a detailed study of the impact 
of DM capture and annihilation on the evolution of Pop. 
Ill stars with a suitable modification of the Geneva stellar 
evolution code 13]. With respect to previous analyses, 
this already allows us to properly take into account the 
stellar structure in the calculation of the capture rate, 
that we compute, following Ref. [3], as 
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Figure 1: Evolutionary tracks of a Pop III 20 Mq star for 
different WIMP densities (labels in units of GeV cm -3 ). We 
have adopted a WIMP model with m x = 100 GeV and of D = 
10" 38 cm 2 . 



Figure 2: Temperature of the core as a function of the DM 
density for the 20 Mq model, at different stages of the core H- 
burning phase. X c denotes the mass fraction of hydrogen at 
the centre ( X c = 0.76 at the beginning of the core H-burning 
phase). WIMP parameters as in Fig. 1. 



where Pi(r) is the mass density profile of a given chem- 
ical element in the interior of the star and Mj refers to 
its atomic mass, while p x , m x and v are respectively the 
WIMP mass and the WIMP density and velocity dis- 
persion at the star position. The velocity of the star 
with respect to an observer, labeled as v*, i s assumed to 
be equal to V, giving therefore r\ = \/3/2. The radial 
escape velocity profile depends on M(r), i.e. the mass 
enclosed within a radius r, v 2 (r) — 2 GM{r )/r 2 dr' . 
The WIMP scattering cross section off nuclei, <t x ,jv is 
constrained by direct detection experiments and for a 
WIMP mass of 100 GeV the current upper limits are 
a s i = 10" 43 cm 2 [H| and a SD = 10" 38 cm 2 respec- 
tively for spin-independent and spin-dependent WIMP 
interactions off a proton. We will adopt these reference 
values throughout the paper, but the capture rate can be 
easily rescaled for other scattering cross sections by using 
Eq. [3] The spin-independent interactions with nucleons 
inside nuclei add up coherently giving an enhancement 
factor A A with respect to the interaction with a single 
nucleon: (J^ N = A A o x . v , where A is the mass number. 
There is no such enhancement for the spin-dependent 
interactions. We consider the contribution to the cap- 
ture rate from WIMP-hydrogen spin dependent interac- 
tions and WIMP-hclium 4 He spin-independent interac- 
tions, neglecting the presence of other elements because 
of their very low abundance. The contribution of He- 
lium is found to be negligible with respect to that from 
hydrogen. 

Once captured, WIMPs get redistributed in the inte- 



rior of the star reaching, in a characteristic time r t h, a 
thermal distribution [17j : 

n ( r ) = n Q e^ with r x = J 3 Tc (3) 

y 27rGp c m x 

with T c and p c referring to the core temperature and 
density. The distribution results quite concentrated to- 
ward the center of the star: e.g. for a 2OM star im- 
mersed in a WIMP density of p x = 10 9 GeV cm~ 3 at 
the beginning of the core H-burning phase we obtain 
r x = 2 x 10 9 cm, a value much lower than the radius 
of the star, i?* = 10 11 cm. This consideration under- 
lines the importance of an accurate spatial resolution in 
the core to properly treat the luminosity produce from 
WIMPs annihilations. We have also checked that re- 
gardless the extremely high concentrations of WIMPs 
obtained at the center of the stars, the gravity due to 
WIMPs is completely negligible. The number of scat- 
tering events needed for DM particles to thermalize with 
the nuclei in the star is of order ps m x /MH, thus an 
upper limit on the thermalization time can be obtained 
as T t h = {m x /MH)/ (crsDnnv) where njj is the average 
density on the star. 

The WIMPs luminosity is simply L x (r) = 
4ir(av)m x c 2 n x (r). For the annihilation cross sec- 
tion times relative velocity (crv), we assume the value 
3 x 10~ 26 cm 2 , as appropriate for a thermal WIMP, but 
note that the total WIMP luminosity at equilibrium 
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Figure 3: Variation of the core H-burning lifetime as a func- 
tion of the WIMP densities for the Pop III 20 and 200 M 
models. WIMP parameters as in Fig. 1. 
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Figure 4: ZAMS positions of 20 and 200 M Pop. Ill stars 
in the g vs. T e // plane for different DM densities (labels in 
units of GeV cm -3 ). Big red circles correspond to the critical 
WIMP density (see text). WIMP parameters as in Fig. 1. 



does not depend on this quantity. After a time 



an equilibrium between capture and annihilation is es- 
tablished, and this incidentally allows to determine the 
normalization constant hq above. 

We have checked that the two transients t x and T t h re- 
main much smaller, during the evolution of the star, than 
the Kclvin-Hclmotz timcscale, tkh and the timescale 
needed for the nuclear reactions to burn an hydrogen 
fraction AX C — 0.002 of the convective core, r„ uc ;: 

_ GM 2 _ q c AX c M» 0.007c 2 

-£L^: J-J Hf, J-J s|s 

where the * labels quantities relative to the the star 
and q c is the core convective mass fraction. This ar- 
gument justifies the assumption of equilibrium between 
capture and annihilation and the use of the radial dis- 
tribution in Eq. [3] We assume here an average WIMP 
velocity v = 10 Km s -1 , the virial velocity in an halo 
of 10 5 - 1O 6 M at z=20. As for the DM density, semi- 
analytic computations of the adiabatic contraction of DM 
halos [12], [l8[ , in agreement with the results extrapolated 
from simulations of first star formation [l9[ , suggest DM 
densities of order 10 12 GeV cm -3 or even higher. 

We have implemented the effects of WIMPs annihila- 
tion in the Geneva stellar evolution code (see Ref. [l3[ 
for details), and followed the evolution of a 20M Q and 



200-Mq stars for different DM densities. We show in 
FigtH the evolutionary tracks for the 20M Q model, and 
show for comparison (black line) the case of a standard 
Pop III star without WIMPs. For DM densities smaller 
than 10 9 GeV cm -3 the evolutionary tracks closely fol- 
low that of a normal star and they are not shown for 
simplicity. The position of the star at the beginning of 
the core H-burning phase (zero-age main sequence, or 
ZAMS) is obtained when, after a short transient, the lu- 
minosity produced at the center of the star equals the 
total luminosity and the star settles down in a station- 
ary regime. For increasing DM densities the WIMPs lu- 
minosity produced at the center overwhelms the lumi- 
nosity from nuclear reactions and makes the star inflate, 
producing therefore a substantial decrease of the effec- 
tive temperature and a moderate decrease of the star 
luminosity at the ZAMS position, with respect to the 
standard scenario. For p x = 10 10 GeV cm -3 , the en- 
ergy produced by WIMPs present in the star at a given 
time, estimated as E x ~ L x tkh, is, at the ZAMS, ~ 0.8 
times the gravitational potential energy of the star, and 
the star therefore starts to contract. In this phase, the 
core temperature, and consequently also the nuclear re- 
actions, increase. When the latter become comparable 
with the WIMPs luminosity, the standard situation is 
recovered and the evolutionary track joins the classical 
tracks of a star without WIMPs. An important differ- 
ence from standard evolution is that in the first phase, the 
nuclear reactions are slowed down and therefore the core 
H-burning lifetime is prolonged. For Dark Matter densi- 
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ties p x < 1.6 10 10 GeV cm -3 , the picture is qualitatively 
the same, and for these models we only show in Fig.Q]the 
first phases of the evolution. In Fig. [2j we show the core 
temperature as a function of the DM density, at different 
stages of the core H-burning phase. At high DM densities 
hydrogen burns at much lower core temperatures than in 
the usual scenario, till a certain mass fraction is reached, 
e.g. X c = 0.3 for p x = 10 10 GeV cm -3 , and the standard 
evolutionary track is joined. For increasing DM densities 
the nuclear reaction rate is more and more delayed till the 
contraction of the star is inhibited, due to the high DM 
energy accumulated, and the evolution is frozen. In Figfl] 
for p x = 2 • 10 10 GeV cm" 3 and p x = 3 • 10 10 GeV cm" 3 
the stars seems to remain indefinitely at the ZAMS po- 
sition. In Fig. [3] we show the core H-burning lifetime as 
a function of the DM density. In the case of a 20M Q 
model, for p < 10 10 Gev cm" 3 the core H-burning phase 
is prolonged by less then 10 % but the delay increases 
rapidly for higher DM densities. Extrapolating the curve 
we determine a critical density, p c = 2.5 • 10 10 Gev cm" 3 , 
beyond which the core H-burning lifetime is longer then 
the age of the Universe. All the calculations have been 
repeated for the 200M Q model and we find that both the 
20Mq and 200Mq stars evolutions are stopped for DM 

densities higher than 5.3 • 10 10 ( 10 _jg C m^ )~ 1 GeV cm ~ 3 - 
We have also verified that the results weakly depend on 
the WIMP mass, e.g. the core H-lifetime is modified by 
a factor 0.2% and 5% respectively for m x — 10 GeV and 
m x = 100 GeV, if p x = 10 10 GeV cm" 3 . 

It is remarkable that under these circumstances, frozen 
Pop III stars can survive until the present epoch, and can 
be searched for as an anomalous stellar population. In 
FigU]we show the effective temperature and gravity ac- 
celeration at the surface of these frozen Pop III stars, 
kept in the H-burning phase, for different DM densi- 



ties. Frozen stars would thus appear much bigger and 
with much lower surface temperatures with respect to 
normal stars with the same mass and metallicity. Our 
results are qualitatively consistent with the preliminary 
estimates in [III QjJ and the analysis in 10]. However, 
for a given DM density, we obtain a somewhat longer 
core H-burning lifetime with respect to [Io| , possibly due 
to their use of an approximated expression for the cap- 
ture rate. We have also followed, for selected models, 
the evolution during the core He-burning phase. During 
this evolutionary stage, the Dark Matter luminosity is 
lower than the nuclear reaction luminosity, therefore the 
impact of DM annihilations is found to be rather weak. 
For the 2OM model and for p x = 1.6 • 10 10 GeV cm" 3 
the He-lifetime is prolonged by a factor 1.2, rather than 
a factor 37 found for the H-burning phase for the same 
DM density. 

In conclusion, we have adapted a stellar evolution code 
to the study the evolution of Pop. Ill stars in presence 
of WIMPs. We have shown that above a critical DM 
density, the annihilation of WIMPs captured by Pop. Ill 
stars can dramatically alter the evolution of these ob- 
jects, and prolong their lifetime beyond the age of the 
Universe. We have determined the properties of these 
'frozen' stars, and determined the observational proper- 
ties that may allow to discriminate these objects from 
ordinary stars. 
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for useful discussions. During the completion of this work 
we became aware of a related work done independently 
by Yoon, Iocco and Akiyama [20]. Their results, obtained 
with an independent stellar evolution code, appear to be 
in good agreement with our own. 
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